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_ ., ' Abstract 

> 
£T) ' It has been pointed out in ref. p] that in the z/MSM (Standard Model extended by three 

right-handed neutrinos with masses smaller than the electroweak scale), there is a corner in 

the parameter space where CP-violating resonant oscillations among the two heaviest right- 

t-h . handed neutrinos continue to operate below the freeze-out temperature of sphaleron transi- 

tions, leading to a lepton asymmetry which is considerably larger than the baryon asymmetry. 

00 

Consequently, the lightest right-handed ("sterile") neutrinos, which may serve as dark mat- 
ter, are generated through an efficient resonant mechanism proposed by Shi and Fuller [2]. 
We re-compute the dark matter relic density and non-equilibrium momentum distribution 
function in this situation with quantum field theoretic methods and, confronting the results 
with existing astrophysical data, derive bounds on the properties of the lightest right-handed 
neutrinos. Our spectra can be used as an input for structure formation simulations in warm 
dark matter cosmologies, for a Lyman-a analysis of the dark matter distribution on small 
scales, and for studying the properties of haloes of dwarf spheroidal galaxies. 
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1. Introduction 

Ever since the experimental discovery of neutrino mass differences, there has been a com- 
pelling case for the existence of right-handed neutrinos in nature. It turns out, however, to be 
difficult to determine the parameters associated with them with any precision. Indeed, given 
that right-handed neutrinos are gauge singlets, their Lagrangian contains explicit (Majorana) 
mass terms in addition to the usual Yukawa interactions. The known mass differences only 
constrain certain combinations of the Yukawa couplings and Majorana masses, so that the 
absolute scale of the Majorana masses cannot be fixed from the existing data. 

Recently, it has been pointed out [3J H] that if the Majorana masses are chosen to be 
significantly smaller than has been the common choice (this corner of the parameter space 
was named z^MSM, for "neutrino Minimal Standard Model"), then it appears possible to 
find an amazingly complete description of the main cosmological mysteries that cannot be 
explained within the Standard Model. Suppose that there are three generations of the right- 
handed neutrinos, like there are of all other fermions in the Standard Model. Then the lightest 
right-handed, or "sterile" neutrinos, with masses in the keV range, might serve as (warm) 
dark matter [SI [2], [6J- I14) 1 ! : the two heavier right-handed neutrinos, with masses in the GeV 
range and almost degenerate with each other, could account simultaneously for baryogenesis 
and the observed active neutrino mass matrix [16] [3] ; while a non-minimal coupling of the 
Higgs field in this theory to the Ricci scalar might explain inflation [T7j. In fact it can be 
argued that the z/MSM could be a good effective field theory all the way up to the Planck 
scale [18] (for a similar argument in a related theory, see ref. |19j). 

On the quantitative level, though, it is non-trivial to realize all of these possibilities within 
the fMSM. Consider the explanation of dark matter by the lightest sterile neutrinos, for 
instance. There are strong experimental constraints from two sides: from the non-observation 
of an X-ray signal generated by the decay of the dark matter neutrinos on one hand ([20]- 
[31] and references therein), and from structure formation simulations on the other [32]-[36j. 
Combining these experimental constraints with the results of theoretical computations of 
thermal dark matter production due to active-sterile mixing (the so-called Dodelson-Widrow 
mechanism) [5j O J9] [TO] [T2] appears in fact to all but exclude the warm dark matter 
scenario [33]— [36], jgoffi 

Such a negative conclusion is premature, however. First of all, the structure formation 
simulations of refs. |34j-|36] assumed the spectrum of the dark matter sterile neutrinos to be 
thermal, with possibly a modest shift of the average momentum towards the infrared, while in 
reality the deviations from the Fermi-Dirac distribution are substantial [101112] . Second, and 
even more importantly, it has been pointed out in ref. [T] that in the framework of the z^MSM 



1 Various cosmological and astrophysical phenomena related to the dark matter sterile neutrinos have been 
discussed in refs. |15| . 

2 Extending the z/MSM by an extra scalar field allows for an additional mechanism for dark matter sterile 
neutrino creation which evades these limits 37, 38, 39 ; another relaxation follows if the reheating temperature 
after inflation is low (in the MeV range) [401I4T1I42] . 



it is possible to generate a large leptonic chemical potential surviving down to temperatures 
of a few hundred MeV. In this situation the results of the theoretical computation change 
dramatically [21 Ej, and it becomes easier to explain dark matter with sterile neutrinos. 

The purpose of the present paper is to elaborate on the latter possibility. More precisely, we 
re-compute the dark matter relic density in this situation with the quantum field theoretic 
methods introduced in refs. [Ill I12j ; analyse uncertainties related to unknown parameters 
and poorly known QCD phenomena; and compare with previous computations in the lit- 
erature. Our main finding is that if lepton asymmetries in the range n Ue /s > 0.8 x 10~ 5 
exist, where n Ue refers to the asymmetry in electron-like neutrinos and s is the total entropy 
density, then the z/MSM can indeed account for the observed dark matter abundance. This 
bound can be consolidated once structure formation simulations have been repeated with the 
non-equilibrium momentum distributions functions ( "spectra" ) of the sterile neutrinos that 
we derive. In any case, asymmetries in the range n Ve /s > 0.8 x 10 -5 may be reachable in 
the so-called Scenario Ha of parameter values of ref. pQ. Note that the largest possible 
asymmetry leading to successful Big Bang Nucleosynthesis corresponds to a chemical poten- 
tial \hl/T\ < 0.07 at T - 1 MeV (95% CL) [Ml Ej, meaning n u Js < 2.5 x 10" 3 in our 
units (cf. appendix A). The maximal asymmetry which can be produced within the z/MSM 
is somewhat smaller, n Ue /s < 0.7x 10 -3 pQ. 

It is appropriate to stress that even though our considerations are naturally viewed as a 
sequel to the lepton asymmetry generated a la ref. [1] , from a practical point of view the origin 
of the lepton asymmetry plays no role in the present analysis, nor do the parameters related to 
the two heaviest right-handed neutrinos. Indeed the only ingredient entering our computation 
is the absolute value of the lepton asymmetry which, as mentioned, we parametrize through 
the ratio n Ue /s. Recalling that the observed baryon asymmetry is ns/s ~ (0.9 — 1.0) x 
10~ 10 [45, 46j, we hence need to assume a boost of some five orders of magnitude in the 
leptonic sector. Besides z/MSM, another possible origin for such an asymmetry could be the 
Affleck-Dine mechanism [37], if it takes place below the electroweak scale and is based on a 
condensate producing many more leptons than quarks. 

Our presentation is organized as follows. In Sec. [2] we generalize the formalism of refs. [11[ 
[T2] to the charge-asymmetric situation. The resulting equation for the sterile neutrino abun- 
dance is integrated numerically in Sec. [31 and the equation for the sterile neutrino spectrum 
in Sec.[U We discuss the astrophysical consequences of these results in Sec.[5l and conclude in 
Sec. El In appendix A we recall the relations of our characterization of the lepton asymmetry, 
through n Ue /s, to a number of other conventions appearing in the literature. A reader only 
interested in the phenomenological consequences of our analysis could start directly from 
Sec. El 



2. Basic formalism 

Our starting point is the Lagrangian 

C = ^N 1 ipN 1 -^ MiiViiVi - F^ZJcirN! - F^^a L L a + £msm , (2.1) 

where N\ are Majorana spinors, and the subscript "1" refers to the lightest right-handed 
neutrino; repeated indices are summed over; M\ is the Majorana mass that we have chosen 
to be real in this basis; L a are the active lepton doublets; F a \ are elements of a complex 
Yukawa matrix; <fr = iT2(j)* is the conjugate Higgs doublet; and a^ = (1— 7s)/2, cir = (l+75)/2 
are chiral projectors. 

To compute the abundance of N\ from first principles, we make a number of basic as- 
sumptions, following refs. |114I12|. First of all, we restrict to temperatures below a few GeV, 
implying that the electroweak symmetry is broken: (</>) ~ (v/y/2,0), where v — 246 GeV 
is the Higgs field vacuum expectation value. Second, we assume that the mixing angles 
9 2 al = \M D \ 2 al /Mf, where \M D \ al = \vF al \/y/2, are very small, 9 al < 1(T 3 . Then it is 
sufficient to restrict to the leading non-trivial order in a Taylor series in 9^. 

The third assumption concerns the flavour structure of the lepton asymmetry. We will only 
consider the case when the asymmetries in all active species (v e , e^, e#, z/^, /j,l, /j,r, v t -, tl, 
tr) are equal. Strictly speaking, this is not satisfied in the z^MSM, since the generation of 
lepton asymmetries takes place when the reactions changing neutrino flavours freeze out pQ. 
We make this assumption in order to keep the discussion as simple as possible, and also 
because it yields the most conservative constraints, leading to the largest sterile neutrino 
abundance and consequently weakening the X-ray bounds. We discuss the general formalism 
applicable in this setting in Sec. 12. 1[ 

At the same time, there is no reservoir replenishing the lepton asymmetry if a part of it is 
converted to right-handed neutrinos: the CP-violating reactions generating the asymmetry 
cease to take place at temperatures above a few GeV pQ. Together with the third assumption 
this completely fixes the time evolution of the lepton asymmetry; this will be demonstrated 
in Sec. E2J 

The fourth and final assumption asserts that the density of the right-handed neutrinos 
produced is below its equilibrium value. This assumption is necessary for the validity of the 
quantum field theoretic formulation of refs. |114 I12j: on the other hand, it may be violated 
in certain parts of the parameter space. In Sec. 12.31 we outline a phenomenological way to 
correct for a possible violation. 

2.1. Results in terras of a generic lepton asymmetry 

Under the assumption of a chemically equilibrated lepton asymmetry among the active 
species, the formal determination of the sterile neutrino production rate proceeds almost 



exactly as in ref . [TT] , with the only difference that the density matrix of the Minimal Stan- 
dard Model (MSM) now takes the form! 



Pmsm = -^msm ex P[ _ P{Hmsm — ^l-^msm)] . 



(2.2) 



Here (3 = 1/T; pl ^ unlike in ref. [IT] ; and Lmsm is the total lepton number operator 
within the MSM, 



Lmsm 



/ d 3 X ^ [laL 70 Ll + Lr 70 Lr + VaL 70 "aL 

a=e,n,r 



(2.3) 



with l e = e,lfj,^ p, l T = t; and ip L = a L i^, ij) R = a R i). 

According to ref. [UJ, the phase space density /i(£, q) of right-handed neutrinos in either 
helicity state s, 

v dJV 1 ( ' ) (t,x,q) 



obeys the equation 



s=l,2 

d 



Ft- H ^ ]flM = R{T ^ h 



(2.4) 



(2.5) 



Here H is the Hubble parameter, H = dlna(£)/d£, and q% are the spatial components of 
the physical momentum q, defined in a local Minkowskian frame. Repeating the analysis of 
ref. [11] with pl / 0, the source term reads 



1 
R(T,q) = Q ^2\M D \ 2 al Tr{$)a L n F (q° - p L )p aa (Q) + n F (q° + p L )p aa (-Q) a R j , 

(2.6) 
where n F (q) = l/[exp(q/T) + l] is the Fermi distribution function; p aa is the spectral function 
related to the propagator of the active neutrino of generation a; and Q is the on-shell four- 
momentum of the right-handed neutrino, i.e. Q 2 = M 2 . Noting that the solution of Eq. (|2.5p 
only depends on q = |q|, it can be written as [12] 



f°°dTMo(T) f T 



h cS (T) 
h e g(T ) 



(2.7) 



where Mq(T) = Mpi^/A5/A7r 3 g e g(T); g e g{T) parametrizes the energy density e as e = 
7r 2 T 4 g e fj(T)/30; and h c g(T) parametrizes the entropy density s as s = 2-K 2 T 3 h e s{T)/A5. 
Moreover, c 2 is the sound speed squared, given by l/c 2 {T) = 3 + Th' cS (T)/h c ff(T). Note 
that in Eq. (|2.7p . the chemical potential pl may be taken to depend on T in any way, to be 
specified later on from physical considerations. 



3 In general, different chemical potentials have to be introduced for different leptonic flavours. 



To derive an expression for p aa , we follow the steps in Sec. 3.2 of ref. [TT], but without 
assuming anything about the symmetry properties of the active neutrino self-energy for the 
moment. The Euclidean propagator (cf. Eq. (3.10) of ref. [TT]) then becomes 



nf Q (Q) = a L - 



+ i$(-Q) 



%Q -i%{-Q) 

O.R = CLL — = = CLR 



(2i 



We have left out the flavour indices from the active neutrino self-energy £ to compactify the 
notation somewhat, and the tildes are a reminder of Euclidean conventions. 

Carrying out the Wick rotation, we can transform Eq. (|2.8j) into a retarded Minkowskian 
propagator (cf. Eq. (3.11) of ref . [TT]): 



n£,(g u ,q)=lCH<Aq)=a L 



+ ?(-Q) 



Q 2 -2Q-S(-Q) + S 2 (-Q) 



an 



(2.9) 



Writing finally S(<7°±i0 + ,q) = Re£(Q)±iIm£(Q) and correspondingly £(— [g°±i0 + ], — q) 
ReS(-Q) =FiImS(-Q), allows to obtain (cf. Eqs. (3.1), (3.4) of ref. [UJ) 



Paa(Q) 



R t r a 



n^(g U + iO + ,q)-n^(g u -iO+,q 



1 

2i L" 

-Si{rQ)[Q - Re 2 (-Q)] - S B (-Q) Im % (-Q) 

a L 7W, — ^\ , no, — ~ a R 



where 



Sr(-Q) = 
Si(-Q) = 



S R (-Q) + Sj(-Q) 

[Q-ReS(-Q)] 2 -[ImS(-Q)] 2 
-2[Q-ReS(-Q)]-ImS(-Q). 



(2.10) 
(2.11) 



(2.12) 
(2.13) 



These expressions can trivially be written also for p aa (—Q), and the results can then be 
inserted into Eq. (|2.6[) . The outcome constitutes a generalization of Eq. (3.12) of ref. |11] . 

To compactify the resulting equations somewhat, we make the following simplifications. We 
note, first of all, that in the imaginary part ImS, the chemical potential changes the thermal 
distribution functions of the on-shell leptons that appear in the intermediate states. Given 
that we are interested in the case pi/T <ti 1.0, however, these changes are not important, 
and will be ignored in the following. Then we can assume that ImS does not get modified 
by pl and, in particular, that ImS(- Q) = ImS(Q), as is the case for pi = 0. 

As far as ReS is concerned, we note that its general structure can be written as 



Re$ aa (Q) = Qa aa (Q) + ib aa (Q) + fic aa (Q) 



(2.14) 



where u = (1, 0). We note that the function a aa (Q) can be ignored, since it is small compared 
with the tree-level term Q . On the other hand the latter structures in Eq. (|2.14p do not appear 
at tree-level, and need to be kept. We have separated two terms: a function b aa (Q) odd in Q, 
appearing already in the charge symmetric situation, as well as a function c aa (Q), defined to 



be even in Q. The function c aa {Q) must be proportional to the leptonic chemical potential 
(or leptonic net number densities), and it is this function which plays an essential role in the 
following. The explicit expression for c aa (Q) can be extracted from ref. [48]; for q° <C mw, 
we can work to first order in an expansion in 1/rriyy, and then the result reads 



3v / 2GF(l + 4sin 2 f9 v 



n. 



(2.15) 



where Gf = g^ / '^\f2mLr is the Fermi constant and, as mentioned, we assumed that all active 



leptonic densities are equal: n v 



n, 



<--L 



n 



t'B 



n 



/'L 



...0 



With these simplifications, we can write 



Q +Re$(Q) 

Q + ReS(Q)] 2 

Q -Re#(-Q) 

Q-ReS(-Q)! 2 



M 1 2 + 2g°(6 + c) + (6 + c) 2 

Q + ji (b - c) , 

Ml + 2q°(b -c) + (b- c) 2 



(2.16) 
(2.17) 
(2.18) 
(2.19) 



where b = b aa (Q), c = c aa (Q). Furthermore, all appearances of ImS can be written in terms 
of the objects 



I Q = Tr 
I u = Tr 



Q a L Im $ {Q)a R \ = 2 Q • Im £(Q) , 
Jlcil Im$ (Q)clr = 2 u ■ ImE(Q) . 



(2.20) 
(2.21) 



Note, in particular, that ImX has a structure analogous to Eq. (|2.14p . with one term propor- 
tional to Q and another to ft, and that consequently even [ImS] 2 can be written in terms 
of the structures in Eqs. ([230]) . §TTi\ . as 



ImS(Q) 



2 _ -p Q + 2q n Q i u -Miii 

4q 2 



(2.22) 



Inserting these simplifications into Eqs. (12. 6p . (J2.11[)ff2.13j) . we finally obtain 
1 



R(T,, 



(27T) 



3 n 



E i m ^i«i 



a=e,/^,r 



where 



x<jn F (g +/x L ) 

fi«(Q) = 
5/(Q) = 



2Sj(Q)[Mf + g°(6 + c)] - S R (Q)I Q 
S 2 R (Q) + S](Q) 

: M 2 + 2g°(6 + c) + (6 + c) 2 + 
: I Q + (6 + c)/ u . 



+ (c->-c,/X£->-^n , (2.23) 



2r2 



I%-2q»I Q I u + MfI, 

4q2 



(2.24) 
(2.25) 



In general, c aQ = \/2Gf[(1 + 2xw 



"',, 



(1 - 2x w ) ^^ nj^ + 2x w E^ % R + 2n„ a + J2p^ a n ^ 



where a;w = sin 



We remark that the production of the dark matter sterile neutrinos, with masses in the 
keV range, takes place at temperatures below a few GeV (cf. Fig. [2] below). In this case, like 
already at ^l = 0, it is numerically a very good approximation to set the term I u to zero, 
whereby Eqs. (|2.24j) . (|2.25j) simplify further. 

Although the formulae given are valid beyond perturbation theory, a practical application 
does make use of approximate perturbative expressions for the functions 6, c and Iq. It is 
important to realise that at the point of a resonance, where some of the "large" terms (Mf + 
2q°b and ±2q°c) cancel against each other, the magnitude of the remainder is determined by 
higher order terms ((6±c) 2 and Iq). A consistent treatment to a certain order in perturbation 
theory would hence require a correspondingly precise (2-loop) determination of the large 
terms 2q°b, 2q°c. At the same time, in a practical application we are not sitting precisely 
at the point of a resonance, but integrate over its contribution, so that for instance a slight 
misplacement of the precise temperature at which the resonance takes place plays little role. 
Consequently, we continue to use 1-loop expressions for the functions b and c throughout^. 
Nevertheless, as we will see, the fact that at the point of the resonance, Im E plays a role also 
in the denominator, will imply that a large ImS can also lead to a decreased abundance, in 
contrast to the case of non-resonant production, where Im £ essentially only plays a role in 
the numerator. For a recent discussion of various resonance-related phenomena, see ref. [49] . 

To be more quantitative about the role of the resonance, we can work out its contribution 
to the production rate semi-analytically. To a good approximation, the resonance is at the 
point where the function 

F{T) = M\ + 2q°(b - c) (2.26) 

vanishes (this comes from the latter term in Eq. (I2.23p . after the insertion of Eq. (J2.24J) ). 
Around this point, the production rate in Eq. (|2.23h can be approximated as 

R ^ ~ nF li)~4 L) jc r |m ^i [m i + q0{b - c W) + ii (2 - 27) 

* "tlw £ Wd\IiM18{HT)), (2.28) 

where we made use of the fact that Iq is very small, in order to identify a representation of 
the Dirac delta-function. Inspecting the expressions for b and c, the function T is positive 
at very low temperatures (because M^ dominates) and at very large temperatures (because 
b dominates), but for sufficiently large n Ve / s and sufficiently small q/T, the term c overtakes 
the others at intermediate temperatures; there are then two zeros of J-, and it turns out to be 
the lower among these that gives the dominant contribution. We denote the corresponding 
temperature by Tr. Eq. (J2.TJ) can now be integrated, to yield 



fit ^~ V" \M |2 l M (T R )n F (q -fi L ) M\ 
/l( °' ^ ~ Jn. r ' Ul n^ s (T R ) (2vrW \r{T R ) 



OL=e,fl,T 



(2.29) 

HTr)=0 



3 For all quantities not specified explicitly in this section, we use the expressions given in ref. |llj . 



2.2. Time evolution of the lepton asymmetry 

The main formula of the previous section, Eq. (12. 23ft . depends on the leptonic chemical 
potential, fii, and on the lepton asymmetry, n Ue , the two of which are related through 
Eq. (|A.3jl . However, the dependence of these quantities on the time (or temperature) has 
been left open. We now need to insert further physical input in order to fix this dependence. 

It is important to realize, first of all, that no reservoir exists for the lepton asymmetry: as 
explained in ref. pQ, the lepton asymmetry was generated by CP-violating processes active 
at temperatures around a few GeV, which subsequently ceased to operate. 

Second, the mass of the lightest right-handed neutrino, Mi, is much below the temperature, 
Mi -C T. Therefore lepton asymmetry violating processes, whose rate is proportional to Mf, 
can to a very good approximation be neglected. In other words, dark matter sterile neutrinos 
and the active leptons can be characterized by a conserved quantity, which we may call the 
total lepton number. In fact, this physics is effectively already built in in Eq. (12. 6\\ . which 
shows that the rate for generating any of the two sterile neutrino states is a sum of two terms, 
with opposite chemical potentials appearing in them, as is appropriate for "particles" and 
"anti-particles" . 

As a result of these two facts, the resonant transitions from active to sterile neutrinos, or 
more precisely the C-odd part of Eq. (|2.6p . cause the original asymmetry to get depleted. 
If the resonance is very effective, the depletion is fast and thereby rapidly terminates the 
resonance phenomenon. 

To be more quantitative, we make the (optimistic) assumption that the flavour and chi- 
rality changing processes within the active generations are fast enough to stay in thermal 
equilibrium. There is then a reservoir of nine spin-1/2 degrees of freedom (three generations, 
each with left-handed neutrinos and both-handed charged leptons) converting to sterile neu- 
trinos. Denoting the two terms in Eq. (|2.6p by R- and i?+, respectively, Eq. (|2.5p can then 
be split and subsequently completed into a closed set of three equations (we also adopt an 
ansatz removing the terms proportional to the Hubble parameter from Eq. (|2.5p ): 



±{9a 3 (t)n Ue (t)} = a 3 (t)Jd 3 q[R + (T,q)-R^T,q)~\, (2.32) 

where the dark matter spectrum f\ is now represented by the sum fi(T,q) = f-(T,q) + 
f+(T, q). The structure of these equations is such that the total lepton charge in a comoving 
volume, 

Ltot = a 3 (t){9n„ e (t) + /d 3 q [f-(t,q) - f+(t,qj\ } , (2.33) 

indeed remains conserved, as must be the case for Mi — ► 0. 



Note that within the approximation of Eq. (|2.28p . the term R + could be omitted from 
Eqs. (|2.30p - (|2.32p . which would simplify the system somewhat. Another practical simplifica- 
tion is to solve the equations in terms of the temperature rather than the time, as we already 
did in Eq. (|27TD . 

A rough estimate for when the depletion has a substantial impact can be obtained as follows. 
If all of the original lepton asymmetry converts to sterile neutrinos, then rii/s > 9n Ue /s, where 
s is the total entropy density, and 

ni(to) = y"d 3 q/i(io,q)- (2.34) 

Therefore the depletion is substantial if n u Js < m/9s. Evaluating the right-hand side of 
this inequality for the case that sterile neutrinos account for all of dark matter, Eq. f|3.3[) . we 
get n u Js < 4.0 x 1(T 4 x keV/9Mi, i.e. 



< 45 ( 10 G ^ ) . (2.35) 



keV 



In other words, for a small initial asymmetry, the depletion has a significant impact at all 
the masses, while for a large initial asymmetry, the effect of the depletion is subdominant 
(because there is more to deplete), unless the mass is small. 

2.3. Back reaction and equilibration 

The derivation of the formulae that our work is based upon, Eqs. (|2.3U|) - (|2.32|) . contains the 
assumption that the particles produced do not thermalize, i.e., that their density remains 
below the equilibrium value at all times. Let us investigate the validity of this assumption. 

It is relatively easy to establish that the total number density of the sterile neutrinos 
produced does remain significantly below the equilibrium value. Indeed, the density of the 
sterile neutrinos is constrained from above by Eq. (|3.3p . and consequently 

^§r < ".OxlO- 4 ^^ - 9.64xl<r>MTo)^ 

n eq (r ) n eq (T ) Mi M l 

keV 
< °- 072 ^> ( 2 - 36 ) 

where we inserted n eq (To) = 3£(3)Tq /2it 2 (cf. Eq. (|A.4p ) as well as h e g(To) < 75 correspond- 
ing to To < 1 GeV [3D]. Thereby the lightest sterile neutrinos appear indeed to be out of 
equilibrium in the whole range Mi/keV > 0.1 that we are interested in. 

It must be realized, however, that the inequality (12.361) is not sufficient to guarantee the 
absence of problems. Indeed, the spectrum of the sterile neutrinos produced is strongly tilted 
towards the infrared 0E]. Given that we are considering fermions, it must then be checked 
that the Pauli exclusion principle is not violated at small momenta. Although an exact 
quantum field theoretic treatment would guarantee this automatically, the assumption of a 



non-thermal result in the derivation of Eq. (|2.23[) means that this consideration now enters 
as an additional ingredient. We refer to the dynamics that prevents an excessive growth of 
the fermionic density at small momenta as "back reaction" . 

Motivated by Boltzmann equations, and recalling our normalization (cf. Eq. (|2.4j) ). we 
expect that the way that back reaction works is to modify the source terms for the distribution 
functions /_,/+ (cf. Eqs. (gSD, (123UD . dMID) by replacing 

nFJq ~ Ml) n F (q° - fi L ) _ n F (q° + \x L ) n F (q° + /j, L ) 

(2vr) 3 ""* (2vr) 3 /_ ' (2vr) 3 ~* (2vr)3 /+ ' [ 6<> 

However, since this recipe would be purely phenomenological at this stage, and since the 
resulting equations are quite difficult to solve numericalljo, we follow a simpler approach 
in the following. Indeed, we first solve Eqs. (j2.30fl (|2.32j) without back reaction, yielding 
the distribution functions which we denote by /_ , f\ . Subsequently, we construct the 
approximants 

n F (q TIJ.L) A0) 

^ y*)* • It , , 

/T - n F ( g o TML ) (0) • V- M > 

This amounts to a rough iterative solution of the structures suggested by Eq. (|2.37|) : guar- 
antees that /=|= never exceed the equilibrium distributions n F (q° =F Ml)/(2vt) 3 ; and, for f T <^ 
n F(q° T A i L)/(27r) 3 , yields the correct result /^ = /^ . 

In order to estimate the practical importance of the back reaction, we have determined our 
main observables (cf. Sees. El H|) for a number of parameter values both from /I and /^. 
We return to the corresponding error estimates in connection with the numerical data. 



3. Sterile neutrino abundance 

The task now is to evaluate the integral in Eq. (|2.7p . with the lepton asymmetry evolved 
according to Eq. (|2.32[) . producing the distribution function fx(to,q), with q = |q|; and 
then to integrate over q, to get the total number density n\ (the order of the integrations 
over T and q can of course be interchanged). We choose to to be the time corresponding 
to To = 1 MeV, below which active neutrinos start to decouple. In order to present the 
numerical results, we start by introducing some further notation. 

First of all, it is conventional to express the mixing angle as sin 2Q a \. For the very small 
Yukawa couplings that we are interested in, it is an excellent approximation to write 



sin 2 2^ = Wl x = 4^i . (3.1) 



6 If the spectra / T do not appear in R^ on the right-hand side, we can integrate over momenta in Eqs. (|2.30[) . 
(|2.31(l . to obtain a coupled set of ordinary differential equations for integrated densities; on the contrary, if 
the spectra do appear in i?=p on the right-hand side, modes with different momenta q couple to each other and 
need to be solved simultaneously, which makes the problem significantly harder. 
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Second, we introduce a total mixing angle as 



sin 2 2fl = J2 49 2 al , (3.2) 



which is the quantity appearing in the X-ray constraints to be discussed below (cf. Fig. H|). 
Now, we can write the total right-handed neutrino density as ni(to) = ^2 a=e UT n a \{t 







where n a \ is the contribution from active flavour a to the dark matter abundance. This 
contribution can conveniently be characterized through the yield parameter Y a \ = n a \/s. 
The corresponding relative energy fraction is Q a i = M\n a \j 'p cr = MiY a i/(p CT /s). Inserting 
p cr /s ss 3.65 x 10 _9 /i 2 GeV from Particle Data Group [35], and noting that VL a \h 2 can amount 
to at most the experimentally known dark-matter density, f2,-t m /i 2 = 0.1143 ± 0.0034 (68% 
CL) [46], we obtain an upper bound on Y a \\ 

keV 
Y al < 4.0 x 10- 4 x — . (3.3) 

Since Y a \ depends monotonously (though non-linearly, because of the depletion discussed in 
Sec. 12. 2|) on sin 2 20 Q ,i, this equation yields an upper bound on the mixing angle. 

Following ref. [12], we concentrate particularly on two flavour structures. In the non- 
resonant case, for a fixed mixing angle, a hierarchy Y e \ > Y^i > Y T \ can be observed [12] . 
because heavier scatterers suppress the production rate. The largest abundance, and the 
most stringent upper bound on sin 2 2#, is then obtained with 

\M D \ el / , \Md\uo_ = \M D \ Tl = "case 1" , (3.4) 

while the smallest relic abundance and the weakest upper bound on sin 2 2# is obtained when 

\M D \ Tl + , \M D \ el = \M D \,a = "case 2" . (3.5) 

As already mentioned, in the resonant case the roles of what leads to the strongest and weakest 
upper bound may get interchanged, because the structure of Eq. (12. 23ft is fairly complicated, 
but it appears that these two cases should still capture the most extreme possibilities. 

To integrate Eqs. (|2.7|) . (|2.32|) in practice, we set the upper limit of the T-integration to 
T max = 4 GeV where R(T, q) is vanishingly small, and the lower limit to To = 1 MeV. We first 
integrate over q, and then evolve n a i, Y a \ and n Ue /s through coupled ordinary differential 
equations in T. This is repeated with several sin 2 2# in order to find the value satisfying the 
constraint in Eq. (|3.3p . Our numerical implementation follows that in ref. [12]. In particular, 
as mentioned above, ImS can be evaluated just as in the case without a lepton asymmetry. 
At the same time, the existence of a narrow resonance does make the integrations over q and 
T much more demanding than in the charge-symmetric case; most importantly, the resolution 
in the g-direction needs to be significantly increased. 

We also remark that at small Mi/keV and large asymmetries, a direct numerical integra- 
tion becomes increasingly difficult, but at the same time Eq. (|2.29p becomes more accurate. 
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Figure 1: Examples of the T-evolution of the lepton asymmetry n Vli j s (cf. Sec. 12. 2 [I . for a fixed 
Mi = 3 keV. Left: a = e. Right: a = t. Note that our results differ even qualitatively from ref. [Hj 
where the asymmetry crosses zero at some temperature. 



However, the asymmetry gets rapidly depleted in this regime, so that in fact Eq. (|2.29j) is a 
good approximation only at the early stages of the resonance. We have found that a work- 
able method is to approximate R as a sum of a C-odd and C-even part; the C-odd part is 
approximated by Eq. (|2.29p . while the C-even part, which dominates for small asymmetries, 
is approximated by the full R from ref. [12], with fiz equal to zero. We have checked that 
the results obtained this way extrapolate, within our resolution, to the "exact" results which 
can be reliably determined at large masses, Mi > 10 keV. 

In Fig. Q] we show examples of the evolution of the lepton asymmetry for various initial 
values. It can be observed that the resonance is quite narrow, and quite effective; in partic- 
ular, for n Uf ,j ' s < 10~ 6 , most of the initial lepton asymmetry is rapidly converted to sterile 
neutrinos, so that the resonance becomes ineffective. Sterile neutrinos are then dominantly 
produced thermally, like in ref. [12], and the mixing angle needs to be large. If the lepton 
asymmetry reservoir were smaller, for instance with three components rather than nine, then 
the depletion would be even more rapid. 

In Fig. [2] we show the resonance temperatures (where existent) for two momenta and 
various asymmetries, as a function of the mass M\. We note that for M\ of a few keV, the 
production peaks at temperatures very close to the QCD crossover. This introduces severe 
hadronic uncertainties to the results, as will be discussed below. 

In Fig. Owe show the upper bound on the mixing angle following from Eq. (13. 3D . for 
n Ue /s = 16.3 x 10~ 6 . This value has been chosen in order to allow for a comparison with 
fig. 1 of ref. [51]. It can be seen that at large masses, M\ > 3 keV, the general order of 
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Figure 2: The resonance temperature corresponding to Eq. (|2.29p . for the modes q/To = 1 and 
q/To = 3, with To = 1 MeV. Left: a = e. Right: a = r. It is seen that, for a given Mi, the resonance 
first affects the smallest values of q/To, and that the resonance extends to larger Mi with increasing 
asymmetry (the asymmetry is indicated in units of 10 6 n I/(! /s on top of the curves). 



magnitude of our result is remarkably close to the result of ref. [51], despite the fact that 
this work was using a simplified kinetic equation and different approximations. On the other 
hand, at small masses, where the depletion discussed in Sec. 12.21 is effective, the results are 
dramatically different. 

In Fig. [3] we have also considered the effects of two sources of hadronic uncertainties: from 
the equation-of-state, which is defined to correspond to a 20% rescaling of the pseudocrit- 
ical temperature of the QCD crossover; and from hadronic scatterings, which is defined to 
correspond to evaluating the hadronic contributions to the vector and axial current spectral 
functions with non-interacting quarks, and an effective number N c = or N c = 3 of colours. 
For justification and more details on this phenomenological but nevertheless conservative 
procedure, we refer to Sec. 5 of ref. [12]. For plotting the dashed and dotted lines in Fig. [31 
we have simultaneously set both uncertainties to their maximal values. It is seen that the 
resulting error depends strongly on parameters, but can be as large as 50%. 

The theoretical upper bound from Eq. (13. 3p is compared with experimental constraints 
(from the non-observation of any X-ray sterile neutrino decay peak in various presumed dark 
matter concentrations) in Fig. HI A more detailed discussion concerning the implications of 
Fig. H] follows in Sec. 

Finally, we note that the plots in this section were produced without taking into account the 
back reaction discussed in Sec. 12.31 i.e., by using the distributions /i . By recomputing Y a \ 
for a number of masses and asymmetries from the corrected distributions f-j- , we find that the 
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Figure 3: The parameter values that, according to our theoretical computation, lead to the correct 
dark matter abundance in the Shi- Fuller scenario [2J; if additional sources are present, sin 28 must 
lie below the curves shown (cf. Eq. (|3.3|) ). For better visibility, the results have been multiplied by 
Mi/keV. The grey region between case 1 (lower solid line on the left, upper solid line in the middle and 
on the right) and case 2 (other solid line) corresponds to different patterns of the active-sterile mixing 
angles, cf. Eqs. (13. 4|l . 1I3.5|) . The dotted and dashed lines correspond to one of these limiting patterns 
with simultaneously the uncertainties from the equation-of-state and from hadronic scatterings set to 
their maximal values. The thick dotted line marked with "Abazajian et al" shows the result in Fig. 1 
of ref. jMj (the case L = 0.003). 



errors are maximal, ~ 25%, for small masses, Mi < 3 keV, and intermediate asymmetries, 
n u e /s ~ 10 — 20 x 10 -6 . For larger masses and other asymmetries, the error from the omission 
of the back reaction is typically below 10%, which we estimate to be well below our other 
systematic uncertainties. 



4. Sterile neutrino spectrum 



We now move from the integrated total sterile neutrino abundance, Eq. ()2.34|) . to the mo- 
mentum distribution function, Eq. (|2.7|) . The physics context where this plays a role is 
that of structure formation, particularly at the smallest scales (Lyman-ct data). The cor- 
responding constraints are considered to be subject to more uncertainties than the X-ray 
bounds, both as far as direct observational issues are concerned, as well as with regard to 
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Figure 4: The central region of Fig. [3j Mi = 0.3 . . . 100.0 keV, compared with regions excluded 
by various X-ray constraints [22J, [25j [30l [31] , coming from XMM- Newton observations of the Large 
Magellanic Cloud (LMC), the Milky Way (MW), and the Andromeda galaxy (M31). SPI marks the 
constraints from 5 years of observations of the Milky Way galactic center by the SPI spectrometer on 
board the Integral observatory. 



dark matter simulations, which have not been carried out with actual non-equilibrium spec- 
tra so far. Nevertheless, adopting a simple recipe for estimating the non-equilibrium effects 
(cf. Eq. (15. ID ), the results of refs. [MIES] can be re-interpreted as the constraints Mi > 11.6 
keV and Mi > 8 keV, respectively (95% CL), at vanishing asymmetry [12]. Very recently 
limits stronger by a factor 2-3 have been reported [36]. We return to how the constraints 
change in the case of a non-zero lepton asymmetry in Sec. We note, however, that the 
most conservative bound, the so-called Tremaine-Gunn bound |52t I53j. is much weaker and 
reads M\ > 0.3 keV |54j . which we have chosen as the lower end of the horizontal axes in 
Figs, m El 

In Fig. [5] we show examples of the spectra, for a relatively small mass Mi = 3 keV (like 
in Fig. [1]), at which point the significant changes caused by the asymmetry can be clearly 
identified. The general pattern to be observed in Fig. [5] is that for a small asymmetry, the 
distribution function is boosted only at very small momenta. Quantities like the average 
momentum {q) s then decrease, as can be seen in Fig. [6j For large asymmetry, the resonance 
affects all q; the total abundance is strongly enhanced with respect to the case without a 
resonance, but the shape of the distribution function is less distorted than at small asymmetry, 
so that the average momentum (q) s returns back towards the value in the non-resonant case. 
Therefore, for any given mass, we can observe a minimal value of (q) s in Fig.[6l (q) s > 0.3(q) a . 
This minimal value is remarkably independent of Mi , but the value of asymmetry at which 
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Figure 5: The distribution function f a i(to,q), for To = 1 MeV and Mi = 3 keV, normalised to the 
massless equilibrium value, feq(to,q) = 2riF(<z)/(27r) 3 . Left: case 1. Right: case 2. These results can 
be compared with refs. [3J \7\: the general feature of strong enhancement at small momenta is the 
same, but our distribution functions show more structure. The case n Vr j s = 16 x 10~ 6 is particularly 
complicated (and sensitive to uncertainties), since the resonance happens to lie just on top of the 
QCD crossover, at T ~ 150 - 200 MeV, cf. Figs. EH3] 

it is reached decreases with increasing Mi. 

Let us stress, however, that the values in Fig. [6] were produced without taking into account 
the back reaction discussed in Sec. 12.3} i.e., by using the distributions /^ . By recomputing 
(q) s /(q)a f° r a number of masses and asymmetries from the approximatively corrected distri- 
butions fzp (shown in Fig. E]), we find that for small masses, M\ < 3 keV, and intermediate 
asymmetries, n Ue /s ~ 10 — 20 x 10~ 6 , the results in Fig. [6] may be too small by up to ~ 25%. 
Thus, for small masses, the minimal average momentum may be better approximated as 
(q) s > 0A(q) a . For larger masses and other asymmetries, the error from the omission of the 
back reaction is typically below 10%, which we estimate to be below our other systematic 
uncertainties. 



5. Astrophysical constraints 

The purpose of the present section is to combine the results of the previous two sections, and 
derive the astrophysical constraints that follow from them. Let us start by briefly recapitu- 
lating, once more, the two different types of considerations that we have carried out. 

First of all, the theoretical computation described in Sec. [3] produces, for a given mass M\ 
and mixing angle sin 2 2#, a definite total abundance of sterile neutrinos. Requiring that this 



16 



1.0 






case 1 






1 


i i i i| i 


1 





' 




-(10 






10 6 n v /s 




£ 0.8 






^^ 


* 

KM 
^ 25o 


" 


II 
H°0.6 

A 

* 0.4 

V 










"- 


A" 

v 0.2 
n n 




, 


1 




" 



10 



10 
M, /keV 



10" 



> 



1.0 



0.8 



case 2 



= 0.6 



A 

* 0.4 
v 



A 
5S< 



0.2- 



0.0 



^— ""- 






"1 

10 6 n /s 

V 
2500 


- 




■7 


' .? 


w 




a 










'n 










•:< 










- 






I.I 


- 



10 



10 

Mj / keV 



10 



Figure 6: The average sterile neutrino momentum, (q) s , normalised to the active neutrino equilibrium 
value, (q) a = 7tt 4 T /180C(3) sw 3.15T . Left: case 1. Right: case 2. 

abundance account for all of the observed energy density in dark matter, leads to the (lepton 
asymmetry dependent) mass-angle relation shown in Fig. HI 

The most direct constraint on the sterile neutrino dark matter scenario comes from com- 
paring these curves with X-ray observations (Fig. |3J). For n Ve / s = 0.0 we are in the allowed 
region only for M\ < 3 keV. Increasing the asymmetry to n Vf ,/s ~ 8 x 10~ 6 opens suddenly 
a whole range of allowed mass values, up to M\ ~ 25 keV. Increasing the asymmetry further 
relaxes the upper bound even more but rather slowly; for instance, if the asymmetry is in- 
creased to n Ue /s ~ 25 x 10 -6 , then we read from Fig. 2] the upper bound M\ < 40 keV, while 
the maximal allowed asymmetry n Ue / s ~ 2500 x 10~ 6 yields the upper bound M\ < 50 keV. 

Another important effect comes from the modification of the sterile neutrino spectrum 
through a lepton asymmetry. As already found in ref. [2], the non-equilibrium spectrum 
of the dark matter sterile neutrinos created in the presence of a lepton asymmetry is very 
different from the thermal one. Some examples are shown in Fig. [5l 

Now, an observation of small scale structures in the Lyman-a data puts an upper bound 
on the free-streaming length and, consequently, on the average velocity of the dark matter 
particles. This converts to a lower bound on the inverse velocity, which, in the absence of an 
actual analysis with non-equilibrium spectra, can be roughly estimated as |32j 



Mi 






> Mo O M x > M ( 



(q)c 



(5.1) 



where (q) a and (q) s are the average momenta of active and sterile neutrinos, respectively, at 
the moment of structure formation, and the value of Mq is Mq ~ 14 keV (95% CL) according 
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to ref. [E! (or M ~ 10 keV at 99.9% CL). According to ref. [35] the bound is somewhat 
weaker, while according to ref. [36] it could be as strong as Mq ~ 28 keV (95% CL). Let us 
start by considering the most conservative bound Mq = 10 keVLj 

The dependence of {q) s /(q)a on Mi and the lepton asymmetry is shown in Fig. [H Quite 
interestingly, this ratio can decrease to 0.3 (or 0.4 at small Mi, cf. end of Sec. H|) for a 
certain range of asymmetries. However, {q) s /{q) a does not decrease further with increasing 
asymmetry, but increases again. Therefore, the lower limit Mq = 10 keV corresponds to 
Mi > 4 keV. 

Combining now the two constraints (from X-rays, Fig. 01 and from structure formation, 
Eq. (|5.ip ). we observe that a solution satisfying both constraints exists for n Ue /s > 8x 10~ 6 . 
The solution corresponds to masses Mi ~ 4 — 25 keV and mixing angles sin 2 2# ~ 8 x 10 -10 — 
2 x 10 -12 . If the lepton asymmetry is increased, larger masses and smaller mixing angles 
become possible. 

Let us then consider the case Mq ~ 28 keV [36]. Using the minimal value (q) s /(q}a ~ 0.3, 
Mq ~ 28 keV corresponds to Mi > 28 x 0.3 ~ 8.4 keV according to Eq. (|5.1j) . Combining 
this with the X-ray constraints in Fig. 0] we see that this can in fact be satisfied with the 
same asymmetry as before, n Ue /s > 8 x 10 -6 . 

It should be stressed, however, that the validity of Eq. (|5.1jl for spectra as extreme as those 
in Fig. [5] remains to be cross-checked. Nevertheless, if true, we can establish the absolute lower 
bound Mi > 4 keV for sterile neutrinos capable of accounting for all of the dark matter in 
the Universe (assuming that the only mechanism for dark matter sterile neutrino production 
is active-sterile mixing). 



6. Conclusions 

The z/MSM, i.e. Minimal Standard Model extended by three right-handed neutrinos with 
masses smaller than the electroweak scale, has a number of parameters not appearing in the 
Standard Model: three Majorana masses and a 3 x 3 complex matrix of Yukawa couplings. 
In ref. pQ, the part of the parameter space associated with the two heaviest right-handed 
neutrinos was explored in detail, and a phenomenologically interesting corner was identified. 
Specifically, it was found that if the mass difference of the heavy (dominantly right-handed) 
neutrino mass eigenstates is much smaller than the known mass differences of the light (domi- 
nantly left-handed) mass eigenstates (Scenario Ha), then it is possible to explain the known 
active neutrino mass differences and mixings, and simultaneously generate and subsequently 
maintain a significant lepton asymmetry in the Universe, without violating constraints related 
to Big Bang Nucleosynthesis at temperatures of about 0.1 MeV. 



T In this work we only consider lower bounds on the mass of the dark matter particle from structure 
formation. However, the problems of missing satellites and cuspy profiles in Cold Dark Matter cosmological 
models, as well as that of the galactic angular momentum, suggest that an upper bound may exist as well [55] . 
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The purpose of the present paper has been to constrain the parameters associated with the 
lightest of the right-handed neutrinos, referred to with the subscript "1". In this case there 
are no constraints from the known active neutrino mass differences and mixings; rather, the 
contribution from the lightest right-handed neutrinos to the see-saw formulae is much below 
0.01 eV. (Consequently the i^MSM excludes the case of degenerate active neutrinos with 
a common mass scale 3> 0.01 eV; in particular, the effective mass for neutrinoless double 
beta decay cannot exceed 0.05 eV |56j.) In contrast, the assertion that all of dark matter 
be made of these neutrinos does allow us to place further constraints on the parameters. 
More precisely, we were led in Sec. [5] to the mass range Mi ~ 4 ... 50 keV. The lower bound 
originates from combining the theoretical analysis of the present paper with observational 
X-ray (Fig. U|) and structure formation (Eq. (|5.ip . Fig. [6|) constraints, whereas the upper 
bound is dictated by the maximal lepton asymmetry allowed by Big Bang Nucleosynthesis. 
The absolute values of the Yukawa couplings of the lightest right-handed neutrinos should be 
in the range 5 x 10~ 15 ...Ax 10~ 13 in this case. 

Of course, these constraints are relaxed if the sterile neutrinos only account for a fraction of 
the dark matter (see, e.g., ref. [57]); if a part of them are produced by some non-equilibrium 
mechanism not related to active-sterile mixing (see, e.g., refs. [37 } 138 ^ 139]): or if the thermal 
history of the Universe is non-standard (see, e.g., refs. [4"0~1 [441 [32] ) . 

It is important to stress, in addition, that the lower bound Mi > 4 keV relies on a naive re- 
interpretation of structure formation simulations which were carried out assuming a thermal 
spectrum of dark matter particles, rather than a proper non-equilibrium shape as given in 
Fig. [5j Hopefully this issue can be put on more solid ground soon. 

Finally, we recall that perhaps the most realistic hope for an experimental detection of dark 
matter sterile neutrinos in the parameter range that we have discussed, would be through the 
discovery of a peak in the diffuse X-ray background from regions where dark matter decays. 
The dominant decay channel is N\ — > v^j and the spectrum should thus peak at the energy 
Mi/2 > 2 keV. As far as laboratory searches are concerned, they are quite difficult due to 
the very small Yukawa couplings of the dark matter sterile neutrinos; however, a possibility 
has been suggested in ref. [58] . 
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Appendix A. Different characterizations of lepton asymmetry 

Several different conventions are used in the literature for characterizing the presence of 
a non-zero lepton asymmetry, and for completeness we specify the relations between them 
here. In analogy with the characterization of the baryon asymmetry through ub/s, the 
conceptually cleanest way is to give the ratio of lepton asymmetry density over the total 
entropy density, because this quantity remains constant as a function of the temperature 
(as long as the Universe remains in thermodynamic equilibrium). However there are many 
leptonic species, so we need to specify which ones to count; in this paper we have assumed 
that the asymmetries are equal in all active leptons (both left-handed and right-handed), and 
choose the electron-like neutrinos (two degrees of freedom, particles and antiparticles) as a 
representative. The corresponding asymmetry, "particles minus antiparticles" , normalised to 
the entropy density, is then denoted by n Ue /s. 

Another way to express the asymmetry is to give the leptonic chemical potential, //£, that 
corresponds to n Ue . This is useful, since chemical potential is the quantity that appears in 
quantum field theoretic computations. In the free limit, 

n Ve = (peLl^eh) = 2i E --5-7 — 2 = / T^W u f( e ~ ^l) ~n F (E + ^ L ) , (A.l) 

where Q{ are the fermionic Matsubara momenta, Qf = (cjf, q) with Wf = 2-7rT(n + ^ ), n € Z; 
Q = (<?Cbq) = Qi + (ifJ'LjO); and E = -\/q 2 + m?. In the massless limit m <§; T, Eq. (jA.lj) 
evaluates to 



Therefore, for /ul -C T, 



^ K ^_ x W , (A.3) 

s 4ir 2 h e g T 

where /i e flF parametrizes the total entropy density through s = 2-7r 2 T 3 /i e fj/45. 

Yet another way to characterize the asymmetry is to use the "intuitive" measure of the 
asymmetry that one is naturally lead to in Boltzmann-equation or density-matrix type treat- 
ments, where the basic degrees of freedom are on-shell particle states: 

A^|^|^ = ^, (A.4) 

where n eq = #^ e + # Pe = 2 / d 3 q/(2vr) 3 /(e g/T + 1) = 3C(3)T 3 /2vr 2 . Then we get 

Or, one can normalize with respect to the photon density, 

L = ^ , (A.6) 

n 7 
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where n 7 = 2/d 3 q/(2vr) 3 /(e 9 / T - 1) = 2C(3)T 3 /vr 2 , which yields 

n Ve 45C(3) 



S 7T 4 /l eff 



x L . (A.7) 



Numerically, choosing to cite the temperature-dependent asymmetries (/j,l/T,A,L) at T 
100 GeV as e.g. in refs. [Hi ETj , we may approximate h e Q ~ 102 for the MSM|f| whereby 



^ « 3.7 x 10- 3 x ^ 
s T 

-3 



(Ai 

T = 100 GeV 



4.1xl0- 3 x A| T = 100GeV (A.9) 

5.4xl0- 3 x L| T = 100GeV . (A.10) 
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